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The charged-current quasi-elastic scattering of muon neutrinos on calcium and 
argon targets is calculated for neutrino energy up to 2.8 GeV. The calculations are 
done within the framework of the relativistic distorted-wave impulse approximation, 
which was earlier successfully applied to describe electron-nucleus data. The model 
is first tested against experimental data for electron scattering off calcium and then 
it is applied to calculate (anti)neutrino cross sections on 40 Ca and 40 Ar. We show 
that reduced exclusive cross sections for neutrino and electron scattering are similar. 
A significant nuclear model dependence of both inclusive and total cross sections 
for energy about 1 GeV was found. From the comparison of the (anti)neutrino 
differential and total cross sections per (proton)neutron, calculated for the carbon, 
oxygen, and argon targets it is evident that the cross sections decrease slowly with 
the mass-number of the target due to nuclear effects. 

PACS numbers: 25.30.-c, 25.30.Bf, 25.30.Pt, 13.15.+g 



I. INTRODUCTION 

The investigation of neutrino oscillations and the precision measurements of neutrino os- 
cillation parameters brought extremely intense neutrino beams. In this situation, statistical 
uncertainties are negligible compared to those systematic uncertainties in the neutrino flux, 
neutrino-nucleus cross sections, and detector effects on both neutrino events selection and 
neutrino energy reconstruction. Thus, the discovery and study of neutrino oscillations have 
renewed interest in neutrino-nucleus interactions. The high intensity neutrino beams allow 
to study these processes with unprecedented details. 

In order to study neutrino oscillation effects in accelerator-based experiments, the neu- 
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trino beams cover the energy range from a few hundred MeV to several GeV. In this en- 
ergy range, the dominant contribution to the neutrino-nucleus cross section comes from 
the charged-current quasielastic (CCQE) reaction and resonance production processes. The 
CCQE is the simplest interaction that represents a two-particle scattering process with single 
final-state proton and it may form a two track events sample. In this process the neutrino 
energy may be estimated using kinematic or calorimetric reconstruction. 

The criteria used to select CCQE events is strongly influenced by both the target material 
and the detector technology. In fully active, fine-grained detectors with good resolution, 
selection techniques that relay on the identification of a single final state proton and lepton 
can be applied. One option for these detector designs is to make their active elements 
as water [1], scintillator [2-5] bars or drift chambers [6] formed in to detector planes. A 
continuous series of these planes forms a region that serves as both primary target and 
tracking detector. 

Another option for fully active and fine-grained detectors is liquid argon time projection 
chamber (LAr TPC) [7, 8] which offers a relatively good particle identification. The high 
spatial resolution and energy measurement down to the MeV scale provides information for 
low and high energy particles. This information allows reduced background for the events 
of interest and potentially improved cross section measurements. The LAr TPC detectors 
are well-suited for long baseline v e appearance physics because of their high efficiency for u e 
'signal' events and low background from events [9]. Therefore, there is a growing interest 
in measuring neutrino cross sections on argon. 

Unfortunately, the cross section data for lepton scattering on argon in relevant energy 
range are rather scarce. There are only experimental data for inclusive 700-MeV electron 
scattering off 40 Ar [10]. On the other hand the structures of 40 Ca and 40 Ar nuclei are almost 
identical and for calcium precise measurements have been performed and are described 
below. 

High resolution exclusive (e, e'p) experiments, where a proton is emitted with a direct 
knockout mechanism on 40 Ca and 48 Ca were carried out at Tokyo [11, 12], Saclay [13], and 
NIKHEF [14-16]. Specific quantum numbers and spectroscopic factors have been assigned 
to the peaks in the observed energy spectrum by studying the missing energy and momentum 
dependence of experimental cross sections. The data analysis of this processes was performed 
within the theoretical framework of the nonrelativistic distorted- wave impulse approximation 
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(DWIA) [17, 18] and relativistic DWIA (RDWIA) [19] in Refs. [14-16, 20-25]. There, 
the RDWIA approach was able to describe with high degree of accuracy the experimental 
shape of the outgoing particle momentum distributions. In order to reproduce experimental 
cross sections, normalization of the bound-state wave functions were fitted to the data and 
identified with the spectroscopic factors. 

Inclusive (e, e') cross sections of the electron scattering on calcium were measured with 
good accuracy in a series of experiments [26-30]. In Refs. [28, 30] the transverse and longi- 
tudinal nuclear response functions were extracted. The comparison of different models with 
the data is described in Refs [31, 32]. 

In this work we compute the exclusive, inclusive and total cross sections for the CCQE 
neutrino scattering from 40 Ca and 40 Ar using the RDWIA model. The calculations of the 
inclusive cross sections are performed with our approach, which includes the final state 
interaction (FSI) effects in the presence of short-range nucleon-nucleon (NN) correlations 
in the ground state [33]. This approach was successfully applied in Refs. [34-37]. First 
we compare our model in describing 40 Ca(e, e'p) 39 K and 40 Ca(e, e') data. Then we apply it 
to the calculation of the CCQE cross sections for the neutrino scattering on 40 Ca and 40 Ar 
nuclei. 

The goals of this work are the following: (a) calculation of the RDWIA CCQE z/ 40 Ar 
cross sections, (b) comparison of the total cross sections, scaled with the number of (pro- 
ton)neutrons in the target for (anti) neutrino scattering on the carbon, oxygen, and argon 
targets, and (c) investigation of nuclear effects on these cross sections. 

The paper is organized as follows. In Sec. II we present briefly the formalism for the 
CCQE scattering process and basic aspects of the model used for the calculation. The results 
are presented and discussed in Sec. III. Our conclusions are summarized in Sec. IV. 

II. FORMALISM OF QUASI-ELASTIC SCATTERING AND RDWIA 

In this section we consider shortly the formalism used to describe electron and neutrino 
quasi-elastic exclusive 



1(h) + A(p A ) -+ l'{k f ) + N(p x ) + B(p B ), 



(1) 
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and inclusive 

l(ki) + A(p A )->lXk f )+X (2) 

scattering off nuclei in the one-photon (W-boson) exchange approximation. Here I labels 
the incident lepton [electron or muon (anti) neutrino], and I' represents the scattered lepton 
(electron or muon), fcj = (£j, fej) and kf = (ef, kf) are the initial and final lepton momenta, 
Pa = {^AiPa)i and Vb = (sb,Pb) are the initial and final target momenta, p x = (e x ,p x ) is 
the ejectile nucleon momentum, q = (ou, q) is the momentum transfer carried by the virtual 
photon (W-boson), and Q 2 = —q 2 = q 2 — u 2 is the photon (W-boson) virtuality. 



A. CCQE neutrino-nucleus cross sections 

In the laboratory frame, the differential cross section for the exclusive electron (cr el ) 
and (anti) neutrino CCQE (a cc ) scattering, in which only a single discrete state or narrow 
resonance of the target is excited, can be written as 

d5ad = JPJ£* £/ « 2 j (el) w nu(el) /O A 

de f dn f dtt x (2tt) 3 Si Q 4 { } 

d 5<jCC _ R \Px\ £ x \kf \ G 2 COS 2 ^c L (cc) W ixv(cc) ^b) 



defdVtfdVt x (27r) 5 £j 
where R is a recoil factor, Qf is the solid angle for the lepton momentum, Q x is the solid 
angle for the ejectile nucleon momentum, a ~ 1/137 is the fine-structure constant, G ~ 
1.16639 xKT 11 MeV" 2 is the Fermi constant, 9 C is the Cabbibo angle (cos 9 C 0.9749), 
L^y is the lepton tensor and and wjf^ are the electromagnetic and weak CC nuclear 

tensors, respectively. The energy e x is the solution to the equation 

e x + e B - m A - u = 0, (4) 

where e B = ^/rn 2 B + p 2 B , p B = q — p x , P x = \/ £ x — m 2 , and rriA, mg, and m are masses of 
the target, recoil nucleus and nucleon, respectively. The missing momentum p m and missing 
energy e m are defined by 

P m = P x ~ Q (5a) 
e m = m + m B - m A - (5b) 

The leptonic tensor is separated into a symmetrical and an anti-symmetrical components 
that are written as in Ref. [33]. For electron scattering at energies considered in this work, a 
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simple and accurate enough effective momentum approximation [23, 38] to include the effects 
of Coulomb distorted wave functions is used. All information about the nuclear structure 
and FSI effects is contained in the electromagnetic and weak CC hadronic tensors, 

(cc) 

and W^v , which are given by the bilinear products of the transition matrix elements of the 
nuclear electromagnetic or CC operator jjf l ^ cc ^ between the initial nucleus state \A) and the 
final state \Bf) as 

Kt l)(cc) = E( s /'^i J i e0(cc) i A )^i J ^ )(cc)t i s /'^)' ( 6 ) 
/ 

where the sum is taken over undetected states. 

The experimental data of the (e, e'p) reaction are usually presented in terms of the reduced 
cross section 

where K el = Rp x e x /(2n) 3 and K cc = Rp x e x /(2n) 5 are phase-space factors for electron 
and neutrino scattering and aw is the corresponding elementary cross section for the lepton 
scattering from the moving free nucleon. The reduced cross section is an interesting quantity 
that can be regarded as the nucleon momentum distribution modified by FSI. It was shown in 
Refs. [33, 35] that reduced cross sections for (anti)neutrino scattering off carbon and oxygen 
are similar to the electron scattering apart from small differences at low beam energy due 
to effects of Coulomb distortion of the incoming electron wave function. 



B. Models 

We describe the lepton-nuclear scattering in the impulse approximation (IA), assuming 
that the incoming lepton interacts with only one nucleon of the target, which is subsequently 
emitted, while the remaining (A-l) nucleons in the target are spectators. When the nuclear 
current is written as a sum of single- nucleon currents, the nuclear matrix element in Eq.(6) 
can write as 

(p,B\J»\A) = J d 3 r exp(it-r)* (_) (p,r)r^(r), (8) 

where T M is the vertex function, t = e B q/W is the recoil-corrected momentum transfer, 
W = a/ (m A + u) 2 — q 2 is the invariant mass, $ and are the relativistic bound-state 
and outgoing wave functions. 
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For electron scattering, we use the CC2 electromagnetic vertex function for a free nu- 
cleon [40] 

r" = 4 d >(QV + ^ v ^f { m\Q\ (9) 

where = ?[7 M , 7 !/ ]/2, Fy 1 ^ and are the Dirac and Pauli nucleon form factors. The 
single-nucleon charged current has V— A structure J^ cc ) = J v + J^. For a free-nucleon 
vertex function r^ cc ) = T v + T A we use the CC2 vector current vertex function 

T v = F V {Q 2 )^ + ^|^F M (Q 2 ) (10) 

and the axial current vertex function 

= F A (Q V75 + Fp(QV75- (11) 

The weak vector form factors F v and F M are related to the corresponding electromagnetic 
form factors Fy 1 ^ and F^ for protons and neutrons by the hypothesis of the conserved vector 
current. We use the approximation of Ref. [39] for the Dirac and Pauli nucleon form factors. 
Because the bound nucleons are the off-shell we employ the de Forest prescription [40] 
and Coulomb gauge for the off-shell vector current vertex T v . The vector-axial Fa and 
pseudoscalar Fp form factors are parametrized using a dipole approximation: 

Fa{q) - {i + QVM\r p{Q) ~ <+Q 2 ' ( } 

where Fa(0) = 1.267, Ma is the axial mass, which controls Q 2 -dependence of Fa(Q 2 ), and 
is the pion mass. 

In RDWIA calculations the independent particle shell model (IPSM) is assumed for 
the nuclear structure. In Eq.(8) the relativistic bound-state wave function for nucleons $ 
are four-spinors that are obtained as the self-consistent (Hartree-Bogolioubov) solutions 
of a Dirac equation, derived within a relativistic mean-field approach, from a Lagrangian 
containing a, uj, and p mesons [41-43]. We use the nucleon bound-state functions calculated 
by the TIMORA code [43] with the normalization factors S(a) relative to full occupancy of 
the IPSM orbitals of 40 Ca and 40 Ar. The source of the reduction of the (e, e'p) spectroscopic 
factors with respect to the mean field values are the short-range and tensor correlations, 
which arise from the characteristics of the bare nucleon-nucleon interaction and long-range 
correlations related to the coupling between single-particle motion and collective surface 
vibrations. 
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For an outgoing nucleon, the simplest choice is to use the plane- wave function \1>, assuming 
that there is no interaction between the ejected nucleon N and the residual nucleus B. In 
this plane-wave impulse approximation (PWIA) the exclusive cross section is factorized into 
the product of the phase-space factor K, elementary off-shell lepton-nucleon scattering cross 
section, and the hole spectral function. Thus, in the PWIA the reduced cross section can 
be interpreted as the momentum distribution of the emitted nucleon when it was inside 
the nucleus. For more realistic description, FSI effects should be taking into account. In 
the RDWIA the distorted wave function \1/ is evaluated as a solution of the Dirac equation 
containing a phenomenological relativistic optical potential. This potential consists of a real 
part, which describes the rescattering of the ejected nucleon and an imaginary part which 
accounts for its absorption into unobserved channels. 

Using the direct Pauli reduction method [24] , the system of two coupled first-order Dirac 
equations can be reduced to a single second-order Schrodinger-like equation for the upper 
component of the Dirac wave function \E'. This equation contains equivalent nonrelativistic 
central and spin-orbit potentials which are functions of the relativistic, energy dependent, 
scalar, and vector optical potentials. We use the LEA program [44] for the numerical 
calculation of the distorted wave functions with the EDAD1 parametrization [45] of the 
relativistic optical potential for calcium. This code was successfully tested against A(e, e'p) 
data for electron scattering off 12 C and 16 [46, 47] and we adopted this program for neutrino 
reaction [33]. 

A complex optical potential with a nonzero imaginary part generally produces an ab- 
sorption of the flux. For the exclusive A(l, I'N) channel this reflects the coupling between 
different open reaction channels. However, for the inclusive reaction, the total flux must be 
conserved. In Refs. [48, 49], it was shown that the inclusive CCQE neutrino cross section of 
the exclusive channel A(l, I'N) is calculated with only the real part of the optical potential 
is almost identical when calculated via the Green's function approach [48], in which the FSI 
effects on inclusive reaction A(l, I'X) are treated by means of a complex potential, and the 
total flux is conserved. We calculate the inclusive and total cross sections with the EDAD1 
relativistic optical potential in which only the real part is included. The inclusive cross sec- 
tions with the FSI effects in the presence of the short-range NN correlations are calculated 
using the method proposed in Ref. [33]. In this approach the contribution of the NN cor- 
related pairs is evaluated in the PWIA. The FSI effects for the high-momentum component 
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are estimated by scaling the PWIA cross section with the A(ef, Clf) function determined in 
Ref. [33]. 

III. RESULTS AND ANALYSIS 
A. Electron scattering 

In this work the IPSM is assumed for the 40 Ca and 40 Ar nuclear structures. The model 
space for 40 Ca(/,/'A) consists of lsi/2, IP3/2, lpi/2, 1^5/2, 2si/ 2 and lds/ 2 nucleon-hole states 
in 39 K and 39 Ca nuclei. The model space for 40 Ar(/, I'N) consists of lsi/ 2 , IP3/2, lpi/2, 1^5/2, 
2sy2 and ld 3 / 2 nucleon-hole states in 39 CI, and lsi/ 2 , lps/2, I-P1/2, 1^5/2, 2sy 2 , 1^3/2 an d 
I/7/2 nucleon-hole states in 39 Ar. All states are regarded as a discrete states even though 
their spreading widths are actually appreciable. 

First we show the performances of the LEA program in describing of experimental data 
for the 40 Ca(e,e'p) reaction measured at NIKHEF in (q,u) constant kinematics [14, 16]. 
The comparison with the experimental reduced cross sections is displayed in Fig.l. We have 
considered cross sections for the removal of the proton from the ld 3 / 2 shell, for transition to 
the l/2 + excited state of the 39 K nucleus at the excitation energy of £^=2.522 MeV, and for 
the transitions to the 5/2 + excited states at E x — 5.258 MeV and E x — 6.328 MeV, obtained 
by knocking out protons from the 2s!/ 2 and ld$/ 2 orbitals, correspondingly. It should be 
noted that positive (negative) values of p m refer to situations where the angle between the 
outgoing proton momentum p m and the incident electron ki is larger (smaller) than the 
angle between q and ki. 

The missing momentum distribution calculated in the RDWIA approach is shown in 
Fig.l with NIKHEF data [14] and provide a good description of the shape of the measured 
distribution. Normalization factors have been applied to reproduce the magnitude of the 
measured reduced cross sections. The factors of 0.68, 0.51, 0.11, and 0.13 for the transition to 
the 3/2+, 1/2+ (.^=2.522 MeV), 5/2+ (£,=5.258 MeV), and 5/2+ (£,=6.328 MeV) states, 
correspondingly are almost identical to those obtained in the data analysis of Refs. [16, 25]. 
Neutrino and antineutrino reduced cross sections of 40 Ca(z/, /i~p) 39 Ca and 40 Ca(z/, ,u + n) 39 K 
reactions also shown are in Fig.l. They were calculated with the same reduced factors as 
electron cross sections. There is an overall good agreement between calculated cross sections, 
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FIG. 1: (Color online) Comparison of the RDWIA calculations for electron (solid line), neutrino 
(dashed line), and antineutrino (dashed-dotted line) reduced cross sections for the removal of 
nucleons from 1^3/2, 2si/2, and W5/2 shells of 40 Ca with NIKHEF data [16]. We show the results 
obtained in (g, u) constant kinematics for electron beam energy £'; )earn =483.2 MeV, the outgoing 
proton kinetic energy T p =100 MeV, and q=450 (MeV/c) 2 . The cross sections are presented as 
functions of missing momentum p m for the transition to the 3/2 + ground state, l/2 + (E x =2.522 
MeV), and 5/2+ (^=5.258, 6.328 MeV) excited states of 39 K and 39 Ca. 
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FIG. 2: (Color online) Reduced cross sections of the Ar(e, e'p) (solid line) and 40 Ca(e, e'p) (dashed 
line) reactions as functions of missing momentum p m for the transitions to the l/2 + (E x =2.522 
MeV) and 5/2 + (E x =6.328 MeV) excited states in 39 C1 and 39 K. For comparison the data for the 
40 Ca(e, e'p) reaction are shown from Ref. [16]. 

but the values of the electron cross sections at the maximum is systematically higher than 
those for (anti)neutrino. This can be attributed to Coulomb distortion upon the incident 
electron wave function. The small difference between neutrino and antineutrino is due to 
difference in the FSI of the proton and neutron with the residual nucleus. 
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TABLE I: Proton and neutron binding energies (E^) and the occupancies for 40 Ca. 



Orbital 




Eh (MeV) 

J - / o v c v / 


11 


S 


lSl/2 


57.4 




64.3 


1 


1P3/2 


36.5 




43.5 


0.95 


lPl/2 


31.6 




38.7 


0.95 


ld 5/2 


15.4 




22.6 


0.80 


2si/ 2 


10.9 




18.1 


0.85 


W3/2 


8.3 




15.6 


0.82 



In Fig. 2 we compare our results for the 40 Ar(e, e'p) reaction for transitions to the l/2 + 
(^=2.522 MeV) and 5/2+ (^=6.328 MeV) excited states of 39 K with the 40 Ca(e, e'p) 
data. Also in this case, the results are multiplied by the same normalization factors (0.51 
and 0.13) as for the 40 Ca(e, e'p) reduced cross sections. In this figure for comparison shown 
are the calculated reduced cross-section for the 40 Ca(e, e'p) reaction. The cross sections for 
the removal protons from the l/2 + and 5/2 + shells of 40 Ca and 40 Ar as functions of p m are 
very similar, but at the maximum the values of the cross sections for 40 Ar are systematically 
lower (less than 12%) than for 40 Ca. 

Mean values of the proton and neutron binding energies and occupancies of shells that 
used in this work are listed in Table I for 40 Ca and in Table II for 40 Ar. The values of the 
proton and neutron binding energies for the Ig^/2, 2si/ 2 , and ld 5 / 2 shells, as well as those 
for neutron for the I/7/2 orbital in 40 Ar were taken from Ref. [50]. For the IP3/2, lpi/2, and 
lsi/2 deeply bound orbitals proton and neutron binding energies were estimated in Refs. 
[51, 52]. In Ref. [16] the 2p 3 / 2 , 1/7/2, 1^3/2; 2si/ 2 , and ld 5 / 2 strengths were obtained as a 
sum of the strengths arising from the discrete transitions and strengths which are observed 
in the continuum at higher excitation energies. In Table I occupancy of the l<i 3 /2, 2s!/ 2 , and 
lo?5/2 orbitals from Ref. [16] are shown. The occupancy of the 1^3/2 and lpi/2 shells were 
estimated from the 40 Ca(e, e'p) data analysis in Ref. [13]. For the lsi/2 shell we assume that 
occupancy is equal to 1. 

Note, that in the IPSM the 20 protons and 20 neutrons of 40 Ca fill the shells up to the 
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TABLE II: Proton and neutron binding 


; energies (Ei) and the 


occupancies for 40 Ar. 


^1 Dllcll 














P 


n 


P 




n 


1 ^ 

J-Si/2 


tW .4 


04. o 


1 




i 
i 


1P3/2 


36.5 


43.5 


0.95 




0.95 


lPl/2 


31.6 


38.7 


0.95 




0.95 


ld 5/2 


15.4 


22.6 


0.80 




0.80 


2si/ 2 


10.9 


18.1 


0.85 




0.85 


ld 3 /2 


8.3 


15.6 


0.85 




0.82 


lf 7/2 




9.87 






0.82 



lc?3/2 orbital whereas the first empty orbital is the I/7/2 orbital. Ground-state correlations 
manifest themselves by a depletion of the orbitals bellow and a filling the orbitals above 
Fermi level, i.e. nucleons from the ld 3 /2 and 2s!/ 2 shells are promoted to the I/7/2 and 
2s 1/2 orbitals. We include the observed l/ 7/2 (^=2.814 MeV) and 2p 3/2 (^=3.019 MeV) 
strength of 0.36 in the 2si/ 2 shell because the observed missing-energy spectrum [16] of 
these orbitals strongly overlap. The rest of the 1/7/2 and 2p 3 / 2 strength were included in 
the ld 3 / 2 shell. In the IPSM the 18 protons (22 neutrons) of 40 Ar fill the shells up to the 
1^3/2 (I/7/2) shell and the occupancies of the shells are not measured. For 40 Ar we assume 
the same occupancies of the orbitals as for 40 Ca because the structures of these nuclei are 
almost identical. 

In our approach the occupancy of the IPSM orbitals of 40 Ca and 40 Ar are approximately 
87% on average. We assume that the missing strength can be attributed to the short-range 
NN correlations, leading to the appearance of the high-momentum and high-energy nucleon 
distribution in the target. We use the general expression [53] for the high- momentum part of 
the spectral function Phm with the parametrization of the momentum distribution for 40 Ca 
taken from Ref. [54]. In our calculations, the spectral function Phm incorporates 13% of the 
IPSM sum-rile limit. The inclusive cross sections with the FSI effects in the presence of the 
short-range AW correlations were calculated using the approach proposed in Ref. [33]. To 
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FIG. 3: (Color online) Inclusive cross section versus energy transfer uj for electron scattering on 
40 Ca. The data are from Ref. [30] for the electron beam energies i? e =408, 545 MeV and scatter- 
ing angle 0=45.5°; -E' e =298, 348 MeV and 0=90°. As shown in the the key, cross sections were 
calculated with the RDWIA, PWIA, and RFGM. 

test our approach, we calculated the inclusive 40 Ca(e, e') cross sections and compared them 
with data from SLAC [26] and Bates [30] experiments. Figures 3 and 4 show measured 
inclusive cross sections as functions of energy transfer as compared to the RDWIA, PWIA, 
and relativistic Fermi gas model (RFGM) calculations with the Fermi momentum pp=2A9 
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FIG. 4: (Color online) Same as Fig. 3, but the data are from Ref. [30] (squares) for E e =628, 681, 
782 MeV, 0=45.5° and from Ref. [26] (triangles) for E e =500 MeV, 6»=60°. 

MeV and nuclear binding energy €{,=33 MeV. These data cover the range of the three- 
momentum transfer (around the peak) from \q\ ^300 MeV/c (beam energy E e =408 MeV 
and scattering angle 9 = 45.5° up to \q\ ^560 MeV/c (£ e =782 MeV, 9 = 45.5°). We note 
that relative to the PWIA results, the generic effect of the FSI reduces the cross section 
value around the QE peak and shifts the peak toward lower values of the energy transfer. 
The peak in the RDWIA calculation occurs at the same energy loss as the data and the 
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value of the calculated cross sections (apart from E e =298 MeV, 9 = 90°) generally agree 
with data within 14%. On the other hand the PWIA and RFGM results systematically 
overestimate the data. 

B. Neutrino scattering 

To study nuclear effects on the Q 2 distribution, we calculated with M^=1.032 GeV the 
inclusive cross sections da/dQ 2 for (anti) neutrino energies e u =0.5, 0.7, 1.2 and 2.5 GeV in 
the RDWIA and RFGM approaches. The results for neutrino and antineutrino scattering 
on calcium and argon are presented in Figs. 5 and 6, respectively, which show da/dQ 2 as 
functions of Q 2 scaled with the number of neutron/proton in the target (cross section per 
neutron/proton). Here, the results for calcium obtained in the RDWIA, are compared with 
cross sections calculated in the RFGM. The cross sections for 40 Ca and 40 Ar are almost 
identical. At the maximum the Fermi gas model results for neutrino (antineutrino) are 
higher than those obtained within the RDWIA. The discrepancy equals to 13% (20%) for 
£^=0.5 GeV and decreases to 7% (12%) for e„=2.5 GeV. The neutrino and antineutrino total 
cross sections for CCQE scattering off 40 Ca and 40 Ar, calculated in the RDWIA and RFGM 
approaches with M^=1.032 GeV are shown in Fig. 7 together with data from Refs. [55-60] . 
Also shown are the total cross sections of the exclusive single-nucleon knock out (u^,fi~p) 
and {u^,ji + n) channels for (anti)neutrino scattering from on-shell nucleons. These reactions 
are source of the CCQE two-particle events in the final states. The cross sections are scaled 
with the number of neutron/proton in the target. 

The ratio between the neutrino cross sections calculated in the RFGM and RDWIA 
decreases with neutrino energy from about 1.26 for e^=0.3 GeV to ssl.16 for e v —l GeV 
and down to ^1.09 for e v =2A GeV. For antineutrino cross sections this ratio is about 2.17 
for £j,=0.3 GeV, 1.3 for e v =\ GeV and 1.16 for e v =2A GeV. The calculated results show 
significant nuclear-model dependence for energy about 2 GeV or less. 

From the experimental data shown in Fig. 7 one can conclude that the CCQE total cross 
sections measured in different experiments can vary by 20-40%. The data have large system- 
atic uncertainties due to the pour knowledge of the background contamination in selected 
events and/or the incoming neutrino flux. Obtaining a reliable estimate of the neutrino flux 
is notoriously difficult and remains a challenge. 
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FIG. 5: (Color online) Inclusive cross section per neutron vs. the four-momentum transfer Q 2 for 
neutrino scattering on 40 Ca and 40 Ar and for the four values of incoming neutrino energy: £,,=0.5, 
0.7, 1.2, and 2.5 GeV. As shown in the key, the cross sections were calculated with the RDWIA 
and RFGM (for calcium) approaches. 

Selection techniques that rely on the identification of a single final state proton (two track 
CCQE events) can improve significantly the purity of the QE sample. Moreover, a simul- 
taneous measurement of both two track and single muon track events [6] allow to constrain 
the systematics associated with the FSI and the neutrino flux. The ratio of the exclusive 
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FIG. 6: (Color online) Same as Fig. 5, but cross section per proton for antineutrino scattering. 

(z/^,// - ^) reaction total cross section to the CCQE total cross section is an attractive quan- 
tity because it is supposed to be rather insensitive to the neutrino flux uncertainty and can 
be especially susceptible to the FSI effects. 

We calculated the R ex = a^ t /a tot ratio, where of^ is the total cross section of the 
(Uft, fiN) reaction for the (anti)neutrino scattering on shell-nucleons on the carbon [35], 
oxygen [33] and argon. The results calculated in the RDWIA are shown in Fig. 8 as func- 
tions of (anti) neutrino energy. This figure clearly shows that the ratio reduces with the 
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FIG. 7: (Color online) Total cross section for CCQE scattering of muon neutrino (upper panel) 
and antineutrino (lower panel) on °Ca and 40 Ar as a function of incoming (anti)neutrino energy. 
The solid and dashed lines are the RDWIA results for calcium and argon, respectively, while the 
dotted line is the RFGM calculation. The dashed-dotted line is the RDWIA result for exclusive 
(v,fj,N) reactions. Data points for different targets are from Refs. [55-60] 
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mass number of the target. The function R ex (e u ) has a maximum in the range e u =0.3-0A 
GeV and decreases slowly with neutrino energy. For carbon (oxygen) (argon) the ratio is 
0.59 (0.55) (0.48) at the maximum and ^0.54 (0.46) (0.43) at e u =2.8 GeV. So, due to the 
FSI the contribution of the exclusive channels reduces slowly with the neutrino energy and 
mass- number of the target. 

Most long base-line neutrino oscillation experiments try to reduce systematic errors in 
the extraction of the oscillation parameters by using near and far detectors. Because these 
detectors are not necessarily of the same target material we estimated the difference between 
the total cross sections per (proton) neutron for the (anti)neutrino CCQE scattering off 
12 C, 16 0, and 40 Ar. The ratio R(e v ) = (cr^)nuci/( cr Zt)nuci (Ar/C ratio) was calculated, 
where the cross sections {(J l tot ) nuc i are scaled with the number of neutron/proton in the 
target. The results obtained in the RFGM and RDWIA are shown in Fig. 9. The ratio 
R(e u = {o? ot )nuci/{vtot)nuci (O / C ratio) calculated in Ref. [35] is also shown for comparison. 

The Fermi gas model predicts almost identical values of af ot and af ot . For neutrino 
(antineutrino) scattering the ratio Ar/C increases from 0.89(0.84) at £,,=0.3 GeV up to 
0.99(0.98) at e u =2.8 GeV. In the RDWIA approach the ratios O/C and Ar/C are lower than 
those calculated in the RFGM. For the neutrino(antineutrino) scattering O/C is 0.88(0.85) 
at £^=0.5 GeV and increases slowly with energy up to 0.93(0.92) at 6^=2.8 GeV. On the 
other hand Ar/C ratio of 0.95(0.75) at e u =0.5 GeV decreases with energy up to 0.88 for the 
neutrino interaction and increases up to 0.88 for the antineutrino scattering at e u =2.8 GeV. 

In the SLAC experiments the inclusive cross sections da/dedQ for electron scattering on 
12 C and 16 [61] as well as on 12 C and 16 [26] were measured in the same kinematical 
conditions. Using these data we calculated the (0/C) e i = (da° /dedVt) nuc i/ \do c /dedQ) nuc i 
and (Ca/C) e i = (da Ca /dedtt) nuc i/(da c /dedil) nud ratios, where the differential cross sections 
(da 1 /dedVL) nucl are scaled with the number of nucleons in the targets. Figures 10 and 11 show 
the measured ratios as functions of energy transfer as compared to the RDWIA calculations 
in the QE peak region. 

There is an agreement between the RDWIA results and the (0/C) e i data within the error 
of the experiments, whereas the uncertainties of the measured (0/C) e i ratios of 5-10% are 
the same order as the predicted effects. On the other hand the calculated (Ca/C) e i ratio 
agree well with data where the observed effect of 15% in the QE peak region is higher than 
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FIG. 8: (Color online) Ratio R ex of the (zv, fiN) reaction total cross section to the CCQE total 
cross section for muon neutrino (upper panel) and antineutrino (lower panel) scattering on 12 C 
(dashed-dotted line), 16 (dashed line), and 40 Ar (solid line) vs incoming (anti)neutrino energy. 
The ratio was calculated in the RDWIA approach. 



experimental errors. Thus the RDWIA model predicts that due to nuclear effects the CCQE 
differential and total cross sections per neutron/proton reduces with the mass number of 
the targets. 

To investigate why the (anti) neutrino CCQE total cross sections per neutron/proton 
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FIG. 9: (Color online) Ratio of the total cross sections per neutron/proton R = O/C (solid and 
dashed lines) [35] and R = Ar/C (dashed-dotted and dotted lines) for CCQE scattering of muon 
neutrino (upper panel) and antineutrino (lower panel) as a function of incoming (anti)neutrino 
energy. As shown in the key, the cross sections were calculated with the RDWIA and RFGM 
approaches. 
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FIG. 10: (Color online) Inclusive cross sections (upper panel) per nucleon and ratio O/C (lower 
panel) vs energy transfer ui for electron scattering on 12 C and 16 0. Data for carbon (filled circles) 
and oxygen (filled triangles) are from Ref. [61] for electron beam energies e e =537, 730 MeV and 
scattering angle # e =37.1 °. The solid line is the RDWIA calculation. 

for 12 C, 16 0, and 40 Ar are dissimilar we estimated the nuclear structure, short-range NN 
correlation, and FSI effects. 

The nuclear structure effects are different in 16 O( 40 Ar) and 12 C cross sections due to 
different nucleon binding energies and momentum distributions for all bound nucleon states 
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FIG. 11: (Color online) Inclusive cross sections (upper panel) per nucleon and ratio Ca/C (lower 
panel) vs. energy transfer u for electron scattering on 12 C and 40 Ca. Data for carbon (filled 
circles) and calcium (filled triangles) are from Ref. [26] for electron beam energies e e =500 MeV 
and scattering angle # e =60°. The solid line is the RDWIA calculation. 

in these nuclei. To estimate these effects we calculated the ratios R str (e u ) = 0/C(Ar/C) 
in the PWIA approach without the NN pair contributions, assuming that the occupancy 
of all nuclear shells S a =l in carbon, oxygen and argon. The nuclear structure effect was 
estimated as A str =l — R str . Then we calculated in the same approach the Rnn(£u) ratios in 
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the presence of the A A correlations in the ground states of nuclei. Note, in our calculations 
the averaged occupancies of the IPSM orbitals of 12 C, 16 0, and 40 Ar are: S c =0.89, S o =0.75, 
and SUr=0.87. The AA correlation effect, i.e. the difference between the cross sections as 
a consequence of different AA pair contributions was calculated as A^^=R str — Rnn- 
The FSI effects are different between the cross sections due to interaction of the outgoing 
nucleons with the different residual nuclei. For example, p + n C for the neutrino or n + n B 
for the antineutrino scattering off the 12 C and p + 15 for the neutrino and n + 15 N for the 
antineutrino scattering off 16 0. These effects were estimated as the difference A FS i=Rnn — 
R, where the ratio R is calculated in the RDWIA approach is shown in Fig. 9 and 10. The 
total nuclear effect is sum A tot =A str + Ajvat + ^fsi — 1 ~~ R- 

The result of the cross section comparison for the (anti) neutrino scattering on carbon 
and oxygen as well as on carbon and argon is displayed in Figs. 12 and 13, correspondingly. 
In these figures the differences A str , Ann, Apsi and A tot are shown as functions of neutrino 
energy. Comparison of the cross sections for 12 C and 16 shows that the nuclear structure 
effects for the neutrino and antineutrino interactions are small and similar. The function 
A str increases slowly with energy up to at £„=2.8 GeV. The AA correlation effects 

decrease with increased (anti) neutrino energy. For the neutrino (antineutrino) scattering 
the function A NN (e u ) reduces from w 8%(7%) at £,,=0.5 GeV up to «4%(4%) at £„=2.8 
GeV. The FSI effects are higher for the antineutrino scattering than those for the neutrino 
interaction. For the neutrino (antineutrino) scattering the function A FS i(e u ) decreases from 
«3%(9%) at £„=0.5 GeV down to «2%(3.5%) at £„=2.8 GeV and the total nuclear effect 
A tot (£„)=0.11(0.17) at £,,=0.5 GeV and reduces down to 0.08(0.09) at £,,=2.8 GeV. 

To conclude, the main sources of the distinction between the carbon and oxygen cross 
sections are different AA pair contributions in the 12 C (11%) and 16 (25%) ground states. 
For the antineutrino cross sections the FSI effects are the same order as AA correlation 
effects at e v > 1 GeV. The precise measurement and accurate calculation of the AA cor- 
relation contributions are important for a reliable estimation of the difference between the 
12 C and 16 O cross sections. 

Comparison of the cross sections for 12 C and 40 Ar shows that the nuclear structure effects 
at £„ >0.5 GeV are identical for neutrino and antineutrino interactions. The function 
A str =0.06 at £„=0.5 GeV and increases slowly with energy up to ^0.07 at £„=2.8 GeV. The 
AA correlation effects of 1-2% are small. The FSI effects are higher for the antineutrino 



25 



0.2 
0.18 
0.1 6 
0.14 
o 0.12 

o o.i 

- 0.08 
0.06 
0.04 
0.02 









0.2 
0.18 
0.1 6 
0.14 
O 0.12 

o o.i 

- 0.08 
0.06 
0.04 
0.02 








0.5 



■ 

— \ \ 


1 1 1 1 1 1 J J 1 1 1 1 1 1 J 1 


1 1 1 1 


Z \ \ 

- \ \ 

- \ \ 


antineutrino 0/C 




1 1 1 1 1 


i i i . ±- r -~ ~~ 1 i i i i 1 i i i 


I I 1 I 



1.5 2 2.5 3 

E„ [GeV] 



FIG. 12: (Color online) Nuclear structure (dashed line), NN correlation (dotted line), FSI (solid 
line), and total nuclear effects (dashed-dotted) (see text) for CCQE scattering of muon neutrino 
(upper panel) and antineutrino (lower panel) on 12 C and 16 vs incoming (anti)neutrino energy. 

scattering. For the neutrino (antineutrino) interaction Apsi is about 0.045(0.13) for £^=0.5 
GeV and reduces with energy down to 0.035(0.06) for £j,=2.8 GeV and the total nuclear 
effect A to t=0.12 (0.18) at e„=0.5 GeV and «0.13(0.13) at e„=2.8 GeV. Thus the nuclear 
structure and FSI effects give the dominant contribution to the difference between the 12 C 
and 40 Ar total cross sections per neutron/proton. 
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FIG. 13: (Color online) Same as Fig. 12 but for (anti)neutrino scattering on 12 C and 40 Ar. 

IV. CONCLUSIONS 

In this work, we study electron and CCQE (anti)neutrino scattering on calcium and argon 
targets in different approaches (PWIA, RDWIA, RFGM). The RDWIA model were widely 
and successfully applied to the analysis of the available electron scattering data over a wide 
range of nuclei. 

First, the reduced cross sections for electron and (anti) neutrino scattering off 40 Ca cal- 
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culated in the RDWIA were tested against 40 Ca(e, e'p) data. We found that the result for 
(anti)neutrino scattering are similar to those for electron scattering (apart from the small 
differences due to the Coulomb correction) and the latter are in good agreement with the 
electron data. Also it was shown that the reduced cross sections for the removal of the 
proton from the shells of 40 Ca and 40 Ar are very similar. 

The inclusive cross sections, calculated in the RDWIA model, which has been modified 
with phenomenological spectroscopic factors and nucleon high-momentum components in 
the target were tested against 40 Ca(e, e') data. The results generally agree within about 
14%. On the other hand, in the QE peak region the RFGM overestimates the value of the 
inclusive cross sections at low momentum transfer and the discrepancy with data reduces 
as this momentum increases. Also, it was shown that the measured and calculated in the 
RDWIA inclusive cross sections per nucleon for electron scattering off 12 C, 16 0, and 40 Ca 
decreases with the mass-number of the target in the QE peak region. 

The CCQE total cross sections for (anti) neutrino scattering on calcium and argon pre- 
dicted by the RFGM are higher than those obtained in the RDWIA and the difference 
decreases with (anti)neutrino energy. The relative contribution of the (z/ M , fiN) channels to 
the CCQE total cross section is lower for heavier nuclei due to the FSI and NN pair effects. 
For (anti)neutrino scattering on carbon, oxygen, and argon we compared the total cross sec- 
tions (scaled with the number of neutrons/protons in the target) and found that the cross 
sections calculated within the RDWIA for 16 O and 40 Ar are lower than those calculated for 
12 C. In the RFGM the cross sections for carbon are practically equal to those for oxygen 
and they are higher by about 2-5% than those for argon at e >0.5 GeV. 

We also studied different sources of the distinction between the 12 C, 16 0, and 40 Ar cross 
sections per neutron/proton. We have found that the difference between the 12 C and 16 
cross sections is mainly due to different NN correlation contributions in the ground states 
of the nuclei and FSI effects. The main sources of the difference between the 12 C and 40 Ar 
cross sections are the nuclear structure and FSI effects. 

Thus the RDWIA approach predicts that the CCQE differential and total cross sections 
per neutron/proton reduces slowly with the mass-number of the target due to the nuclear 
effects. Although the model and theoretical ingredients adopted in the calculations con- 
tain approximations, our results can serve as a useful reference for long base-line neutrino 
oscillation experiments 
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